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Abstract 

Number-state expansions are derived for the squeezed states produced by four-wave 
mixing (modulation interaction and phase conjugation) in fibers. These expansions 
are valid for arbitrary pump-induced coupling and dispersion-induced mismatch 
coefficients. To illustrate their use, formulas are derived for the associated field- 
quadrature and photon-number variances and correlations. 
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1. Introduction 



Parametric devices based on four-wave mixing (FWM) in fibers can generate photon 
pairs for quantum communication experiments [U, Three different types of FWM are 
illustrated in Fig. 1. Modulation interaction (MI) is the degenerate process in which two 
photons from the same pump are destroyed, and signal and idler (sideband) photons are 
created (2tt p — > tt s + 7Tj, where ?Tj represents a photon with frequency Uj). Inverse MI is the 
degenerate process in which two photons from different pumps are destroyed and two signal 
photons are created (n p + n q — > 2n s ). Phase conjugation (PC) is the nondegenerate process 
in which two different pump photons are destroyed and two different sideband photons are 



created (71V 



7T 



7T. S 



7Tj). These processes are reviewed in [31 HI [6] . 
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Figure 1: Frequency diagram for (a) modulation interaction, (b) inverse modulation interaction, 
(c) outer-band phase conjugation and (d) inner-band phase conjugation in a fiber. Long arrows 
denote strong pumps (p and q), whereas short arrows denote weak sidebands (s and i). 

The evolution of an optical system is governed by the (spatial) Schrodinger equation 



iH\ 



where z is distance, d z = d/dz, is the state vector, and the Hamiltonian H depends 
on the creation and destruction operators of the interacting modes (a] and cij, respectively, 
where t denotes a hermitian conjugate). The solution of Eq. ([T]) can be written in the 
input-output form 

Mz)) = U(z)M0)), (2) 
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where the evolution operator U(z) = exp(iHz) is unitary. Notice that U(—z) = TP(z). In 
the Schrodinger picture, the mode operators are constants. 

In the Heisenberg picture, one defines the mode operators aj(z) = W(z)aj(0)U(z), which 
evolve according to the (spatial) Heisenberg equations 

d z cij = i[a,j, H], (3) 

where [ , ] denotes a commutator. In the small-signal (undepleted-pump) regime, the Hamil- 
tonian depends quadratically on the sideband operators, so Eqs. (J3J) are linear in these 
operators. Hence, their solutions can be written in the input-output form 

= T,k\Pjk(z)a k (0) + u jk (z)al(0)}, (4) 

where fijk{z) and Vjk( z ) are transfer functions. 

Any measurable quantity associated with mode j can be written as the expectation value 
of some function of the mode operator and its conjugate. Common examples are the field 
quadrature and photon number. Let ( ) denote an expectation value and F be any function 
that has a Taylor expansion. Then 

(F( aj )) = (iP(z)\F[ aj (0)Mz)) = M0)|F[ fli C0]|tK0)). (5) 

In the Schrodinger picture [first part of Eq. (j3J)], the state vector evolves and the mode 
operator is constant, whereas in the Heisenberg picture [second part of Eq. (jSJ)] the state 
vector is constant and the mode operator evolves. Expectation values of output quantities 
(which involve one or more operators) can be calculated using either picture. 

To model photon-generation experiments, one needs to determine the probabilities of 
measuring different numbers of photons. In this context, the Schrodinger picture is prefer- 
able. In this report, number-state expansions are derived for the squeezed states produced 
by (inverse) MI and PC in fibers. These expansions generalize the standard results [3 [8], 
which do not include the effects of fiber dispersion. As written, they apply to parametric 
processes driven by continuous-wave pumps. However, by defining suitable superposition 
(Schmidt) modes [10] , one can also apply them to parametric processes driven by pulsed 
pumps. 
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2. One- mode squeezed state 

The inverse MI in a fiber is governed by the Hamiltonian 

H = 5ala s + [ 1 (al) 2 + 1 *a 2 s }/2, (6) 

where a s is the destruction operator of the signal mode, 5 is the mismatch coefficient and 7 
is the coupling coefficient. Formulas for these coefficients, which involve the fiber dispersion 
and nonlinearity coefficients, and the the pump amplitudes, are stated in [H [6] . 
By combining Eqs. ([3]) and ([6]), one obtains the evolution equation 

d z a s = i5a s + ryaj. (7) 

The solution of this equation can be written in the input-output form 

a 8 (z) = f i(z)a a (0) + u(z)al(0), (8) 

where the transfer functions 

/jl(z) = cos(fcz) + iSsin(kz)/k, (9) 

v{z) = i r ysm(kz)/k (10) 

and the inverse-MI wavenumber k = (5 2 — I7I 2 ) 1 / 2 . Equations (El) and ( fTUl) are based on 
the assumption that k is real (and the MI is stable). If k = in is imaginary (and the MI is 
unstable), cos(A;2;) is replaced by cos\i(k,z) and sin(kz)/k is replaced by smh(Kz) / k. Notice 
that fi(—z) = fi*(z) and v{—z) = —u(z). The transfer functions also satisfy the auxiliary 
equation 

|/z| 2 - |if = 1. (11) 

By definition, the one-mode squeezed state produced by inverse MI is U(z)\0), where 
U(z) = exp(iHz) and |0) is the one-mode vacuum state. One can facilitate the calculation 
of the output state by rewriting U in normally-ordered form. To do this, one writes 

H = -iK + + 25^3 + 7 *iT_ - 5/2, (12) 

where the operators 

K+ = (al) 2 /2, K_ = a 2 s /2, K 3 = (aja s + a s a\)/A. (13) 
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These operators satisfy the angular-momentum-like commutation relations [K+, K-] = —2K^ 
and [7^3, -K"±] = ±K±. By using a standard operator-ordering theorem, which is proved in 
the Appendix, one finds that 

exp(iHz) = exp(y + K + ) exp(7 3 if 3 ) exp(7__^_), (14) 

where the auxiliary functions 

y+(z) = iysm(kz)/[kcos(kz) — iSsin(kz)), (15) 
y~(z) = iy* sm(kz)/[k cos(kz) — iSsm(kz)], (16) 
73(2) = — 2 log [cos (/cz) — id sin(kz) / k) (17) 



and the (inconsequential) phase factor exp(— i5z/2) was omitted. It follows from Eq. (114p . 
and the identities K± = and K\ = K3, that 7+(— z) = 7* (z), 7_(— z) = y+(z) and 
73 (—z) = 73(2). The auxiliary functions satisfy these conditions. By comparing Eqs. ([15]) - 
([T7|) to Eqs. (jUJ) and ([TO]) , one obtains the compact formulas 

-y+ = v/tS, y- = -v*/v*, 7 3 = -21og(^). (18) 
Hence, if the input is the vacuum state, the output is the squeezed state 



(^) 1/2 „=oV/ 2«n! 

where the basis vectors \2n) = (aj) 2n |0)/[(2n)!] 1 / 2 . Notice that each eigenstate contains an 
even number of photons. 

For the special case in which 5 = (maximal exponential growth), 

/jl = cosh(\y\z) , v = iy sinh(\y\z)/\y\, uf/j,* — iytanh.(\y\z)/\y\ (20) 

and Eq. (TT9]) reduces to the standard result [3 Ej. (To verify this statement, use the 
substitution iyz = —se.) For the complementary case in which 5 = \y\ (transitional linear 
growth) , 

fi = l+i\y\z, v = iyz, v/n* — iyz/(l — i\y\z) (21) 
and Eq. ([TO]) reduces to the result of [TTJ. (Use the substitution \y\z = z.) 
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It is customary (and easy) to calculate the moments of a) and a using the Heisenberg 
picture. However, I will calculate the lower-order moments using the Schrodinger picture, 
to check Eq. (fl9l) and illustrate its use. The zeroth-order moment is (^l^). Let c n be the 
coefficient of |2n) in Eq. (fT9l and Pm = \c n \ 2 be the probability of a 2n-photon state. Then 
I A* I -f 2n. = x 2n (2n)\ / 4 n (n\) 2 , where x = \vjp\- It follows from Eq. fjTT]) and the identity 

Six) ~ 1 - T X2n(2n)l (22) 
[ } ~ (l-x 2 )V2 -2- o 4«(n!) 2 1 J 

that = 1: The state vector (fT9|) is normalized. 

The field quadrature 

g s = (at e ^+a s e-^)/2 1/2 , (23) 

where 0; is the local-oscillator phase, and the quadrature deviation Sq s = q s — (q s ). The 
first-order moment 

(q.) = 0, (24) 

because a^ip) and a\ip) contain only odd-number states, whereas (ip\ contains only even- 
number states. The expectation values (means) of all the odd moments are zero, for the 
same reason. To calculate the quadrature variance (Sq 2 ) (which equals (q 2 )), one needs to 
calculate the inner products of (if)\ and 



^ 2 ty) = (fi*MJ2 2 ^n\2n), (25) 



v a „, 

n=l 



ala 



n=l 
oo 



J22nc n \2n), (26) 



a, a 



3 "'gl 



£(2n + l)c n |2n>, (27) 



n=0 



a 



2| 



(^///*)£(2n + l)c n |2n). (28) 



n=0 



It follows from Eq. (J22D that 



;r 2 



V 2nP n = xd x S(x) = — ^-ttj. (29) 

By using this result and Eq. ( fill to evaluate the inner products, one finds that ((aj) 2 ) = 
/iV*, (aja s ) = |z/| 2 , (a s aj) = |/i| 2 and (a 2 ) = \w. Hence, the quadrature variance 

(Sq 2 ) = (|/i| 2 + 2\fiv\ cos6+ M 2 )/2, (30) 



where the phase difference 9 = M + <p v — 20/. The quadrature variance attains its maximum 
+ |z/|) 2 /2 when 9 = and its minimum — \u\) 2 /2 when 9 = n. In the stable regime 
\v\ 2 is bounded by 7 2 //c 2 , whereas in the unstable regime it is unbounded [Eq. fllOB ]. The 
quadrature is squeezed in both regimes. Equation fl30l) is consistent with the results of [6|fT2]. 
which were obtained using the Heisenberg picture. For the special case in which 5 = 0, it 
reduces to the standard result [7J [8] . 

Now define the photon-number operator n s = a\a s and the number deviation 6n s = 
n s ~ ( n s)- It also follows from Eq. (|22|) that 

\li\(n™) = (xd x ) m S(x). (31) 

By using this result and Eq. ( TTTT) . one finds that 

(n s ) = |z/| 2 , (n 2 s ) = W\ 2 (2\fi\ 2 + M 2 ), (Sn 2 s ) = 2\/2u\ 2 . (32) 

Equations (1321) are consistent with the results of P, [12], which were obtained using the 
Heisenberg picture. For the special case in which 5 = 0, they reduce to the standard results 

EE]. 

3. Two-mode squeezed state 

MI and PC in a fiber are governed by the Hamiltonian 

H = 5(ala s + a\a,j) + 7aja| + 7*a s a i; (33) 

where aj is the destruction operator of mode j (s or i). Formulas for the mismatch and 
coupling coefficients are stated in [3l [5]. By combining Eqs. ([3]) and (|33|) . one obtains the 
evolution equations 

d z a s = i5a s + iya\, (34) 
d z a\ = —iy*a s — iSa\. (35) 

The solutions of these equations can be written in the input-output form 

a s (z) = (i(z)a a (0) +v(z)a\(0), (36) 
a\{z) = u*(z)a s (0)+fi*(z)4(0), (37) 



where the transfer functions were defined in Eqs. ([9]) and (flOl) . and the MI (PC) wavenumber 
k = (5 2 - | 7 | 2 ) 1/2 . 

By definition, the two-mode squeezed state produced by MI (PC) is U(z)\0, 0), where 
|0, 0) is the two-mode vacuum state. One can calculate this state by writing 

H = -fK + + 25K 3 + 7*iiT_ - 5, (38) 

where the operators 

K + = ajaj, K__ = a s cii, K 3 = (aja, + a s a\)/2. (39) 

These operators also satisfy the commutation relations [K + ,K^\ = —2K 3 and [K 3 ,K±] = 
±K±. By using the aforementioned operator-ordering theorem, one can rewrite U in the 
form of Eq. flHj), where the auxiliary functions were defined in Eqs. (jl5j) — (II7j) and the 
(inconsequential) phase factor exp(— iSz) was omitted. Hence, if the input is the vacuum 
state, the output is the squeezed state 

M = ^E(£)V»>- (40) 

Notice that each eigenstate contains an equal number of signal and idler photons. For the 
special cases in which 6 = and S = \ Eq. (JiOj) reduces to the standard result [3 [8] and 
the result of [TTJ, respectively. 

Let c n be the coefficient of \n,n) in Eq. (jlDl) and let P n = \c\\ be a probability. Then 
|/i| 2 P n = y n , where y = |z///i| 2 . It follows from the identity 

1 oo 

s(y) = — = Y.y n (4i) 
1 -y 

that = 1- By combining Eqs. (j40p and (141 p . one can show that 

\»\ 2 (n7 j nT k ) = (yd y r^S(y), (42) 

where rij = a]aj and k ^ j. 

The quadrature and photon-number operators of the signal and idler, and their devia- 
tions, are defined in the same way as the signal operators were defined in Sec. 2. Both 
quadratures 

<?i> = 0, (43) 
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because contains states with equal numbers of signal and idler photons, whereas a,j\ifj) 
and dj\ip) contain states with unequal numbers of signal and idler photons: The photon 
numbers are unbalanced. Most of the operator moments vanish: Only powers of 
a^a\ and aja k are nonzero. To calculate the quadrature variances {5qf) (which equal (gj)) 
and correlation (5qj5qk) (which equals (qjqk)), one needs to calculate the inner products of 
(ib\ and 



a)a\\ 



a)aj\ 



aj a]\ 



= Yl n Cn\n,n), (44) 

n=l 

oo 

= ^nc n \n,n), (45) 
n=i 

oo 

= J2( n + 1 ) c n\n,n), (46) 

n=0 

oo 

ajdklip) = (v//i*)J2(n + l)c n \n,n). (47) 

n=0 

It follows from Eq. (I4"2"j) that (ajaj) = /-tV*, (a]cij) = \v\ 2 , (djdj) = | /x | 2 and (djdk) = V>v. 
By combing these results, one finds that 

H) = (i/f + M 2 )A (48) 
(Sq s 5qi) = \fjLv\ cos 0, (49) 

where the phase difference = <fi^ + <fi u — 20;. Neither of the output modes is squeezed by 
itself. (The quadrature variances are phase independent.) Instead, squeezing is manifested 
as a quadrature correlation, which strengthens with distance. It also follows from Eq. (1421) 
that 

(rij) = \u\ 2 , (n)) = M 2 (|/i| 2 + M 2 ) = (rijUk). (50) 
In turn, it follows from Eqs. (1501) that 

(5n 2 ) = |/iz/| 2 = (Srijdnk)- (51) 

Equations ([4"8"]) - (j51~]) are consistent with the results of [3 [12], which were obtained using the 
Heisenberg picture. For the special case in which 5 = 0, they reduce to the standard results 
[H Further analysis shows that ((rij — rik) m ) = 0: The sideband photon-numbers are 
perfectly correlated, as implied by Eq. (|40l) . 
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4. Decomposition of a two-mode squeezed state 

It was stated in Sec. 1 that one can relate multiple-mode transformations to one-mode 
transformations by defining suitable superposition modes [10] . This statement applies to 
the two-mode transformation discussed in Sec. 3. Define the sum and difference modes 

a± = {a s ±a l )/2 1 ' 2 . (52) 

Then, by making these substitutions in Eq. ( 1331) . one obtains the alternative Hamiltonian 

H = Sala+ + [ 7 (4) 2 + 7*a+]/2 

+ 5ala_ - [7(al) 2 + 7*a 2 _]/2, (53) 

in which the + and — terms are separate. H + is identical to the one-mode Hamiltonian 
flS]), whereas in H_ the coupling coefficient —7 has the opposite sign. Hence, the two-mode 
squeezed state (JID]) is the direct product of two one- mode states of the form (fT9l) . where 
H± = /i and v± = ±u. One can also demonstrate this equivalence directly, by writing 

/i~ V/iV n\ 

= ■i£m"M)l^g|o + ,o_) 

1 » " / zy \" (at) 2fc (-l)»- fc (aL) 2 ("- fc ) 

2-k\{n-k)\ 1 + ' ^ 

1 ^^/^y-»/V (4) a (ol) a , n n N 
= ?SSvJ (y) 2*mn\ |0+ '°- ) " (54) 

The last of Eqs. (|54p has the required form. 

5. Summary 

Parametric devices based on modulation interaction (MI) and phase conjugation (PC) 
in fibers can generate photon pairs for quantum communication experiments. In this re- 
port, number-state expansions were derived for the one- mode squeezed state produced by 
inverse MI [Eq. (j!9p ]. and the two- mode squeezed states produced by MI and PC [Eq. (1401) ]. 
These expansions are valid for arbitrary pump-induced coupling and dispersion-induced mis- 
match coefficients. Hence, they apply to a variety of polarization-dependent parametric pro- 
cesses driven by continuous- wave pumps in strongly-birefringent, randomly-birefringent and 
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rapidly-spun fibers. They also apply to the Schmidt modes that participate in parametric 
processes driven by pulsed pumps. To illustrate their use, formulas were derived for the 
associated field-quadrature and photon- number variances and correlations [Eqs. (1301) . (132jl . 
flU, flUD and (JSH)]. 
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Appendix: Operator-ordering theorem 

The main results of this report, Eqs. ffl9|) and fj40|) . were obtained by the use of an 
operator-ordering theorem (00T). Although such theorems are common in the quantum- 
optics literature [7j Ej, they are not common in the optical-communications literature. Con- 
sequently, in this appendix the OOT ( TH|) will be proved from first principles. 

The proof of this OOT relies on the Baker-Campbell-Hausdorff (BCH) lemma 

oo 

exp(a)6exp(— a) = Y^fa, b] n /n\, (55) 

n=0 

where a and b are operators, and the nth-order commutator [a, b] n is defined recursively: 
[a, b]o = b, [a, b]i = [a, b] and [a, b] n = [a, [a, 6] n _i]. There are two ways to prove this lemma. 
The first (direct) way is to expand both sides of Eq. ( |55l) in Taylor series, and equate the 
coefficients of a n [TT]. The second (elegant) way is to define the function 

F(x) = exp(xa)6exp(— xa). (56) 

It follows from Eq. (J56J) that F'(x) = aF - Fa = [a,F] and F"(x) = a[a, F] - [a,F}a = 
[a, [a, F}}, where F' = dF/dx. By extending this sequence, and using the fact that F(0) = b, 
one finds that 

oo 

F{x) = 52[a,b] n x n /n\ (57) 

ra=0 

The BCH lemma is Eq. (1571) . with x = 1. 

Equation (I14p provides a normally-ordered formula for the Schrodinger evolution-operator 
exp(iHz), where if is a Hamiltonian and z is distance. In this report 

H = jK + + 25K 3 + j*K_, (58) 

where 5 is real, K\. = K T and K\ = K 3 . The ii'-operators satisfy the commutation relations 
[K + ,K-] = —2Kz and [K3, K±] = ±K±. (Formulas for these operators were stated in Sees. 
2 and 3.) Because one can multiply K + and K_ by conjugate phase factors without changing 
the commutation relations, one can simplify the derivation of the OOT by assuming that 7 
is real. Define the function 

G(z) = exp[i( 7 K + + 25^3 + lKJ)z\. (59) 
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Because the i^-operators form a closed set under commutation, one can rewrite Eq. (|59|) in 
the normally-ordered form 

G(z) = exp\p(z)K+] exp[q{z)K 3 ) exp[r(z)KJ\, (60) 

where p, q and r are functions of z (to be determined). It follows from Eq. (1591) that 

G' = i{iK + + 26K 3 + iK„)G, (61) 

where G' = dG/dz. Likewise, it follows from Eq. (IBTfl) that 

G' = {p'K + + q'e pK+ K 3 e- pK+ + r'e pK+ e qK3 K_e- qK3 e- pK+ )G. (62) 

By using lemma fl55l) and the aforementioned commutation relations, one finds that 

e pK+ K 3 e- pK+ = K 3 -pK + , (63) 
e qK3 K_e- qKa = K„e~ q , (64) 
e pK+ K_e- pK+ = K_-2pK 3 +p 2 K + . (65) 



By using these results to simplify Eq. (1621) . and equating the coefficients of K + , K 3 and K_ 
in Eqs. fIBTl) and fl62l) . one obtains the differential equations 

p' — pq' + p 2 {r' e~ q ) = ry, (66) 
q'-2p(r'e~ q ) = 2i5, (67) 
r'e" 9 = i 7 , (68) 

respectively. Equations (|661) - (I6SI) are to be solved, subject to the boundary (initial) condi- 
tions p(0) = 0, q(0) = and r(0) = 0. 

By combining Eqs. fl66l) - fl68|) . one obtains the individual equation 

p' = z( 7 - ( 5 2 / 7 ) + i 7 (p + 5/ 7 ) 2 . (69) 
This equation has the implicit solution 

tan -1 [(tp + 5)/k] — tan -1 [S/k] = ikz, (70) 
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where the parameter k = (7 2 — <5 2 ) 1//2 . By inverting Eq. ([TO]) , one obtains the explicit solution 

p(z) = i^y smh(nz) / [k cosh(Kz) — i8 sinh(Kz)]. (71) 

It is easy to verify that 

q(z) = — 2 log [cosh (kz) — i5 smh(Kz) / k], (72) 
r(z) = i / jsmh(Kz)/[Kcosh.(Kz) — idsmh^Kz)] (73) 

are the solutions of Eqs. fl67j) and fl68l) . respectively. To allow for complex 7, one replaces 7 
by 7* in Eq. ( !73l) and 7 2 by I7I 2 in the formula for k. These results are consistent with the 
formulas for j± and 73 [Eqs. (|T5l)-(|T7jl]. 
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